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Q \ Abstract 

. Efficient selection of emission line galaxies at z ^ 1 by photometric information in wide field surveys 

(-H ' is one of the keys for future spectroscopic surveys to constrain dark energy using the baryon acoustic 

Q^, oscillation (BAO) signature. Here we estimate the Ha and [Oil] line luminosity functions of galaxies at 

z = 0.5-1.7 using a novel approach where multi-wavelength imaging data is used to jointly estimate both 
photometric redshifts and star-formation rates. These photometric estimates of line luminosities at high- 
redshift use the large data sets of the Subaru Deep Field and Subaru XMM-Newton Deep Field (covering 
^1 ~ 1 deg^) and are calibrated with the spectroscopic data of the local Sloan Digital Sky Survey galaxies. 

The derived luminosity functions (especially Ha) are in reasonable agreement with the past estimates 
based on spectroscopic or narrow-band-filter surveys. This dataset is useful for examining the photometric 
selection of target galaxies for BAO surveys because of the large cosmological volume covered and the large 
' number of galaxies with detailed photometric information. We use the sample to derive the photometric 

\^ , and physical properties of emission line galaxies to assist planning for future spectroscopic BAO surveys. 

■ We also show some examples of photometric selection procedures which can efficiently select these emission 

I line galaxies. 

. Key words: cosmology: observations - galaxies: statistics - techniques: photometric 

o ■ 

1. Introduction (2006) for a theoretical review.) A key question is whether 

• • : the dark energy behaves exactly like A, or something more 

. ^ ' The precise measurements of the cosmological param- exotic. These cases can be distinguished obscrvationally 

, eters in the past decade have shown that we live in a by higher precision measurements of the expansion history 

' Universe whose expansion is accelerating (Riess et al. of the Universe. 
5^ , 1998; Perlmutter et al. 1999; Spergel et al. 2007; Komatsu One route to this expansion history is by using the 

et al. 2008). This has led to the paradigm of the ACDM baryon acoustic oscillation (BAO) method which has at- 

universe whose acceleration is caused by a non-zero classi- tracted particular attention in recent years as a promising 

cal 'cosmological constant' (A), equivalent in effect to the new probe of dark energy (Blake & Glazebrook 2003; Seo 

energy density of the vacuum in relativistic quantum me- & Eisenstein 2003; Glazebrook & Blake 2005). The BAO 

chanics. Although this simple model is almost perfectly scale in the Universe serves as a 'standard ruler' whose 

consistent with the cosmological data obtained so far, the length is determined by the sound horizon at the time 

value of A is unexplained and in fact 10^^^ times smaller of recombination and which imprints itself in the clus- 

than plausible physical values. These problems have moti- tering of matter as revealed by the cosmic microwave 

vate the generalisation of the accelerating component as a background (CMB) at early times and the distribution 

'dark energy'. The dark energy may be energy density of a of galaxies at later times. Because the BAO scale (~ 150 

still unknown physical field, or modification of the theory comoving Mpc) is set by straight-forward linear physics 

of gravity and spacetime on the cosmological scale may in the early Universe it is relatively robust against as- 

be required. To reveal the nature of dark energy has the trophysical uncertainties and can be calibrated by CMB 

potential to revolutionize our understanding of the uni- observations (Eisenstein & White 2004). The BAO sig- 

verse and fundamental physics. (See, e.g., Copeland et al. natures have already been seen in the distribution of low- 
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redshift galaxies by the Two Degree Field Galaxy Redshift 
Survey (2dFGRS) (Cole et al. 2005) and Sloan Digital 
Sky Survey (SDSS) (Eiscnstcin ct al. 2005; Tegmark et al. 
2006; Percival et al. 2007; Okumura et al. 2008). 

A clear next step is to detect BAO at higher redshifts, 
which would yield constraints complementary to those of 
low redshift surveys. A number of such BAO surveys are 
planned, including both spectroscopic and photometric 
surveys. In this paper we concentrate on spectroscopic 
BAO surveys at z 1 , which are optimal surveys because 
apparent BAO scales at z ~ 1 are the most sensitive to 
dark energy. However BAO surveys need to cover mini- 
mum volumes of order 1 Gpc^ in order to be effective 
(Blake & Glazebrook 2003) which implies survey areas of 
at least 500-1000 deg^. Since the largest fields of view 
on optical telescopes are of order 1 degree then for sur- 
veys to be feasible exposure times must be kept short (be- 
low one hour). Because of this star forming galaxies with 
emission lines are of particular interest, since we can de- 
termine the redshifts with a shorter exposure time. The 
emission lines useful for such surveys are likely to be the 
[Oll]A3727 forbidden line doublet in optical spectroscopy, 
and Ha A6563 in near-infrared spectroscopy, both of which 
are well-known as representative star formation indicators 
(Kcwlcy ct al. 2004; Moustakas et al. 2006). 

Important issues for such BAO surveys are (1) whether 
there are a sufficient number of emission line galaxies that 
can be detected with a reasonable exposure time, and (2) 
how to select them from photometric survey data as the 
targets of spectroscopic observations. These issues can 
be examined by analysing the existing estimates of the 
luminosity functions of emission line galaxies at z ^ 1, 
based on slitlcss spectroscopic surveys (Yan ct al. 1999; 
Hopkins et al. 2000 in Ha), photometry-selected spec- 
troscopic surveys (Tresse et al. 2002 in Ha; Hogg et al. 
1998; Teplitz et al. 2003; Zhu et al. 2008 in [Oil]), pho- 
tometric surveys using narrow-band filters (Geach et al. 
2008 in Ha, Ly et al. 2007; Takahashi et al. 2007 in 
[Oil]), and Fabry-Perot spectroscopic surveys (Hippelein 
et al. 2003 in [Oil]). However, the statistics of these 
studies (especially Ha) are rather limited, because of the 
small survey area in spectroscopic surveys (^ 60 arcmin^ 
for Ha) and small redshift range in narrow-band surveys 
(Az ^ 0.04). In general BAO surveys are targeting bright 
and low number-density galaxies compared with those tar- 
geted in general studies of galaxy evolution, consequently 
limited statistics may lead to significant uncertainties due 
to small statistics and cosmic variances. Recently Zhu 
et al. (2008), covering a large area (~14,000 galaxies in 
2.45 deg~^), has appeared but this is the only large [Oil] 
survey. 

Here we study these issues based on wide-field broad- 
band photometric data of the Subaru Deep Field (SDF) 
(Kashikawa et al. 2004) and the Subaru XMM-Newton 
Deep Field (SXDF) (Furusawa et al. 2008), with estimat- 
ing redshifts and emission line luminosities from star for- 
mation rate based on the photometric redshift calcula- 
tions. Although the redshift and emission line luminosity 
estimates are less reliable than those of spectroscopic or 



narrow-band surveys, the statistics is greatly improved 
by the SDF and SXDF data, which are the deepest data 
among the surveys covering ^ 1 deg^. Furthermore, deep 
photometric data in a variety of band filters allow us to ex- 
amine the photometric properties of emission line galaxies 
in detail, to infer the required depth for imaging surveys 
and efficient procedures to select targets for spectroscopic 
observations. 

We will present our results as far as possible in a general 
way so that they are useful for any BAO survey planning. 
To complement this we also examine the cases of specific 
BAO survey proposals. Though there are several BAO 
surveys planned at z ^ 1, we examine the two proposed 
BAO surveys as examples. One is the FastSound project, 
which is a proposed BAO survey using the Fiber Multi- 
Object Spectrograph (FMOS), which is a new multi- fiber 
NIR spectrograph for the Subaru Telescope. FMOS has 
had engineering first light in May 2008 (Iwamuro et al. 
2008) and is now undergoing further test observations. A 
BAO survey using Ha emission line galaxies at z ^ 1 has 
been proposed (Glazebrook et al. 2004a). The other is 
a BAO survey based on [O II] emission line galaxies us- 
ing the Wide-field Fiber-fed Multi Object Spectrograph 
(WFMOS) (Bassett et al. 2005), which is a proposed next- 
generation optical spectrograph for the Subaru Telescope 
having up to several thousand fibers. 

This paper is organized as follows. In §2 we describe the 
photometric data sets used in this work and the method- 
ology of photometric redshift calculations. In §3, we es- 
timate the emission line luminosities with calibrations by 
the SDSS spectroscopic galaxy sample. We then calculate 
the Ha and [O II] luminosity functions and compare with 
those in previous studies, in §4. In §5 we discuss impli- 
cations for future spectroscopic BAO surveys. Finally we 
show some examples of color selection of target galaxies 
for spectroscopic BAO galaxies in §6. We summarize our 
results in §7. Throughout this paper, we use the stan- 
dard cosmological parameters of Hq = 70 km s~^ Mpc~^, 
= 0.3, = 0.7, fib = 0.045, as = 0.8 in the ACDM 
universe (Spergel et al. 2007). All magnitudes are given 
in the AB magnitude system. 

2. The Sample and Photometric Redshifts 

2.1. The SDF and SXDF Galaxies 

The SDF is located at R.A. = 13'^ 24"" 38^9, Decl. = 
+27° 29' 25". 9 (J2000), where the Galactic extinction is 
E{B -V)^ 0.017 (Schlegel et al. 1998). All magnitudes 
presented below have been corrected for the Galactic ex- 
tinction. We use the 0.114 deg^ region where both the 
Subaru/Suprime-Cam data in BVRi'z' bands (Kashikawa 
et al. 2004) and UKIRT/WFCAM data (Hayashi et al. 
2007; Motohara et al. 2009, in preparation) in JK bands 
are available. The limiting magnitudes of the optical and 
NIR data are shown in Table 1. The seeing size of optical 
images is 0.98 arcsec FWHM (1 pixel — 0.2 arcscc). We 
use 19,494 galaxies detected both in i' and K bands. The 
requirement for K-hand photometry is because of the need 
for accurate photometric redshifts and stellar contamina- 



No. ] Photometric Ha and [O II] LF of SDF and SXDF Galaxies: Implications for Future BAG Surveys 3 



Table 1. Basic Information of SDF and SXDF 



Field 


Area [deg"^] 


Number 






Limiting 


magnitudes* 


[mag" 












B 


V 


Re 


z' J 


K 


3.6/im 


4.5/im 


SDF 


0.114 


17408 


28.5 


27.7 


27.8 27.4 


26.6 24.0 


24.3 






SXDF 


0.732 


76193 

n J. 


28.1 


27.8 


27.6 27.6 


26.6 24.5 


24.0 


23.1 


22.4 



* The limiting magnitudes are Scr at 2"apcrturc for B, V, Rc, and z' bands, and 3cr at 3". 8 aperture for 3.6 and 
4.5/im bands. For the J and K bands, they are 3cr at 2"aperture (SDF) and 3(j at 2". 3 aperture (SXDF). 



tion removal, however we will show below that this is not 
a significant source of bias. 

The SXDF is centered at R.A. = 2^ 18'" 00^0 and 
Decl. = -5° 00' 00".0 (J2000), where the Galactic 
extinction is E{B — V) = 0.020. We use the photomet- 
ric data in BVRi'z' bands (Furusawa et al. 2008) taken 
by Subaru/Suprime-Cam, in JK bands taken by UKIDSS 
survey (Warren et al. 2007), and in 3.6 and 4.5 /im bands 
taken by SWIRE survey (Lonsdale et al. 2004), in the 
area of 0.732 deg^ where the data in all of these bands 
are available. The details of the catalog will be reported 
elsewhere (Akiyama et al. 2009, in preparation). The lim- 
iting magnitudes are given in Table 1, and the seeing size 
of optical images is about 0.8 arcscc FWHM. As in the 
SDF we require A"-band detection, and consequently we 
use 82,323 galaxies detected both in i' and K bands. 

Stellar objects are removed as follows. For bright ob- 
jects with i' < 21.0, we used the SExtractor (Bertin & 
Arnouts 1996) stellarity and objects are removed when 
the i?-band GLASS STAR index is larger than 0.5. For 
objects fainter than i' = 21.0, we use the Bz'K color- 
color diagram, and objects are removed if {z' ~ K) < 
0.5{B - z') - 0.8 and (z' - K) < 0.3(B - z') - 0.4. (All 
galaxies selected by the detection in i'- and _ftr-bands are 
detected also in B and z' bands.) After the removal of 
stars by these procedures, there remain 17,408 and 76,193 
galaxies in the SDF and SXDF, respectively. 

2.2. Photometric Redshift Estimate 

We use the publicly available code hyper z (Bolzonella et 
al. 2000) for our photometric redshift estimations. Wc use 
7 template models of spectral energy distribution (SED) 
having exponentially decaying star formation history with 
the characteristic time scales of r = 0.1, 0.3, 0.7, 1, 3, 
5, and 15 Gyrs. In addition to this, we added a SED 
template with a constant star formation rate (SFR), and 
hence we have 8 SED templates in total. We use the stel- 
lar spectrum library of Bruzual & Chariot (2003) with the 
solar metallicity {Z ~ Zq ~ 0.02) and the initial mass func- 
tion of Salpeter (1955) in the mass range of 0.1-100 Mq. 
We apply the extinction law of Calzetti et al. (2000) for 
dust extinction ranging Ay = 0.0-2.0 at an interval of 0.2. 
The stellar mass is estimated taking into account instanta- 
neous recycling with a typical value of the returned mass 
fraction from evolved stars R ~ 0.25 (e.g., Nagashima & 
Yoshii 2004), as M* = (1 - i?)il/sFR^int, where AfgFR-int 
is the stellar mass calculated by the integration of SFR 
through the star formation history of a galaxy. It is well 
known that the Salpeter IMF overestimates the number of 
low-mass stars compared to more realistic IMFs, We use 




Fig. 1. Comparison of spectroscopic and photometric red- 
shifts of SXDF galaxies. The solid and dashed lines show 
the lines of Zphot = Zspec and l^phot - Zspcc 1/(1 -I- ^spcc) = 0.2, 
respectively. 

the Salpeter IMF as it is convenient to use with hyper z, 
however we divide our final masses by 1.82 to convert them 
to the Baldry & Glazebrook IMF (Baldry & Glazebrook 
2003) which is more realistic but has a particularly simple 
scaling with respect to the Salpeter IMF (Glazebrook et 
al. 2004b). 

The photometric redshifts (zphot) calculated in this way 
are compared with the spectroscopic redshifts (zspec) of 
550 X-ray undetected galaxies in SXDF having observed 
spectra (Figure 1). As in Ilbcrt et al. (2006), we define the 
redshift accuracy as the median of 1.48| Az|/(1 -I- z), and 
outlier fraction as the fraction of those having |Az|/(l + 
z) > 0.2. The redshift accuracy and outlier fraction in 
Zphot < 2.0 and Zgpec < 2.0 are 0.030 and 6%. This level 
of agreement is sufficient for the purposes of this work 
for most of the galaxies. It can be seen that galaxies 
at Zgpoc ~ 0.3 tend to be estimated as Zphot 3.0. be- 
cause of confusion between the 4000A and Lyman breaks. 
However, since we are considering objects in the range 
0.5 < z < 1.7, this confusion will not affect our results. 
The photometric redshift distributions of SDF and SXDF 
galaxies in the different thresholds of A'-magnitude are 
given in Figure 2. 

In the following sections, we will mainly show the result 
of SXDF galaxies, because the photometric redshift cal- 
culations in SXDF are expected to be more reliable than 
those in SDF due to the advantage of having mid-infrared 
data. 
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Fig. 2. The photometric redshift distributions of SDF and 
SXDF galaxies brighter than the X- magnitude of 22 (top), 
23 (middle), and 24 mags (bottom). The number of SDF 
galaxies multiplied by 2 is shown in each panel. 



3. Line Luminosity Calibration by SDSS 
Galaxies 

Wc estimate emission line luminosities from SFRs of the 
best-fit SED templates in the photometric redshift calcu- 
lations. Our primary goal is not to measure the SFRs 
themselves but rather use them as an intermediate step 
towards the final line luminosities. Accordingly rather 
than use standard conversion factors from line luminos- 
ity to SFR (e.g., Kennicutt 1998), we calculate our own 
conversion factors as part of our calibration process. This 
ensures we empirically obtain the best estimate of the re- 
lation between line luminosities and our photometrically 
estimated SFR, including systematic errors such as dust 
and mctallicity effects. The calibration is done by us- 
ing the large spectroscopic sample of local galaxies of the 
SDSS. 

3.1. The SDSS Sample 

We use the u'g'r'i'z' photometric catalog publicly re- 
leased from the SDSS studies at Max-Planck-Institute for 
Astrophysics / Johns Hopkins University^, which is based 
on the fourth SDSS data release (Adelman-McCarthy et 
al. 2006). We select a sample with 14.5 < Tpotro < 17.77 
and 0.023 < z < 0.399 in order to detect Ha and [Oil] 
lines in the wavelength coverage of SDSS, where rpotro is 
the Petrosian magnitude of r-band. We use the sample 
of star-forming galaxies constructed by Brinchmann et al. 
(2004), and the SDSS data set used in this work includes 
75,757 objects, whose mean redshift is 0.081. 

The SDSS photometric zero points are close to the AB 
convention but they arc not reproduced exactly. We con- 
vert the SDSS photometry system into the AB system 
using the AB correction of Eisenstein et al. (2006). The 
SDSS spectroscopic observations cannot measure the total 
line emissions from galaxies, because of the limited size of 
the fiber diameter 3". We correct this aperture effect by 
the relation between Petrosian and fiber magnitudes as 
described in Hopkins et al. (2003). 

3.2. Calibration 

We use the hyperz on the SDSS u'g'r'i'z' catalog for 
our photometric redshift estimations, as in the case of 
SDF and SXDF. In this work, we are interested in es- 
timating observable line luminosities, rather than the 
true extinction-corrected line luminosity corresponding to 
SFR. Therefore wc compare the SFR that is estimated by 
hyperz but reduced by the magnitudes corresponding to 
extinction (Ahq and Ajonj), with the observed line lumi- 
nosities (LHQ,obs and i[oii],obs)- The amount of extinc- 
tion is calculated from the Ay estimates by hyperz and 
the adopted extinction law. 

The correlation between the two quantities of the SDSS 
galaxies is shown in Fig. 3, for Ha and [Oil] lines. It 
should be noted that we have corrected SFR estimates by 
an additional factor of [dL(zspcc) / dLizphot)]'^ , to allow for 
the difference between the true spectroscopic redshift and 
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Fig. 3. The extincted SFRs estimated by photometric red- 
shift calculations versus H« and [O ll] line luminosities mea- 
sured by spectroscopic data of SDSS galaxies. The solid lines 
are the best-fit linear relation. The scatters in the case of Ha 
and [On] are ±0.24dcx and ±0.35dex, respectively. 

that returned by hyper z. For low-redshift samples such as 
SDSS a small rcdshift difference Az does not significantly 
change the observed SED shape (as AAz is much smaller 
than the band width) and hence the hyper z fit. However it 
has a much stronger effect on derived luminosities which is 
removed by this correction. We cross-checked this method 
of SFR estimates using a sub-sample of a few hundred 
SDSS galaxies where we manually adjusted hyperz to re- 
turn the SED best fit at Zgpoc- We found the SFRs from 
the two approaches consistent to within a factor of two. 
We use the former approach as it is much more efficient 
for large samples with the existing hyperz codebase. 

We see a tight correlation between SFR and Ha lumi- 
nosities, and we obtain the best-fit conversion relation as: 

, / ina.obs \ , / SFRphot \ „ . , 
V[erg s i]y \[Mq yr 

-f-40.85. 

The offset -1-40.85 is different from that of Kennicutt 
(1998) by 0.25 dex, which is not unreasonable consider- 
ing the uncertainties in stellar spectrum libraries and ef- 
ficiency of ionization photon production. 

Next we consider the [O II] emission. As can be seen 
in the right panel of Fig. 3, the correlation with SFR 
is not as good as that of Ha emission, which is almost 
certainly due to the effect of metallicity and extinction 
(Kewley et al. 2004; Moustakas et al. 2006). Hence, follow- 
ing Moustakas et al., we calibrate the conversion law from 
SFR to [O II] luminosities as a function of absolute i?-band 
magnitude (M^ phot) in order to remove the dominant lu- 
minosity dependence of these effects. Figure 4 shows the 
conversion factor, ?/[oii]7 versus absolute i?-band magni- 
tude, phot- Here the conversion factor is defined by 
the following equation: 




-22 -21 -20 -19 -18 



Fig. 4. The conversion factor »?[oii] from SFR to [On] 
line luminosity as a function of absolute B-band magnitude, 
Mb 

phot- The small dots are the SDSS galaxies. The open 
circles and errors show the mean and standard deviation of 
the SDSS galaxy distribution in each A/g phot bin. 




Fig. 5. Comparison of spectroscopic and photometric [On] 
luminosities of SXDF galaxies with available spectra. The 
solid lines show the lines of log (Lfon] ,phot/-f-[Oii] , spec) = 
and ±1.0. The scatter is approximately ±0.5dex. 
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Wc derived 7][oii] from calculating best fit values using 
SDSS data in every M^ phot bins, which are shown by 
open circles in Fig. 4. The best-fit formula of the conver- 
sion factor is given as follows: 

log[??[oii](^/B,phot)] =43.24 + 0.13MB 

One possible criticism of this calibration is that it is 
based on local galaxies and hence might not applicable to 
high-z galaxies (for example due to evolution in metallic- 
ity effects). Wc test this issue by using 88 SXDF galaxies 
with [O II] luminosities estimated from the emission line 
equivalent widths, and Figure 5 shows the comparison be- 
tween photometric and spectroscopic [O II] luminosity es- 
timates. Here we again correct the photometric line lu- 
minosity estimates by the factor of [c?L(zspoc)/dL(2:phot)]^ 
as for the SDSS sample, to remove the error coming from 
the distance estimate. The agreement between the two 
is reasonable though the scatter is significant (~0.5dcx). 
We also note also that Moustakas et al. showed that their 
approach for [O II] to SFR conversion worked reasonably 
well for galaxies at z ~ 1. 

4. Line Luminosity Functions 

J^.l. Ha Luminosity Function 

We construct the Ha luminosity function per unit log- 
arithmic luminosity interval in a given redshift interval 
{zl,Z2)^ by simply dividing the observed number of galax- 
ies by the comoving volume corresponding to the redshift 
interval. We do not use the more sophisticated 1/Vmax 
method (e.g., Yan et al. 1999) to estimate LFs, because we 
cannot determine a clear limiting flux of Ha emission for 
galaxies originally selected by broad band fluxes. Instead, 
we discuss the effect of limiting magnitudes on the LF 
estimates below. 

The derived Ha luminosity functions from SDF data 
(filled triangles) and SXDF data (filled circles) are shown 
in Figure 6, for three different redshift ranges of 0.5 < z < 
1.0 (top), 1.0 < z < 1.4 (middle), and 1.4 < z < 1.7 (bot- 
tom) . We also plot the Ha luminosity functions measured 
by previous studies, after correcting the cosmological pa- 
rameter into those used in this work. The redshift ranges 
of the previous studies are slightly different from ours (see 
figure caption), and we plot the previous data in the panel 
of similar redshift ranges in this figure. 

To examine the effect of limiting magnitudes on the 
LF results, the A"-magnitude distributions of galaxies in 
several bins of Xhq are plotted in Figure 7. Since the K- 
band observation is relatively shallower than the i' band, 
we can examine the effect of the limiting magnitudes by 
comparing these distributions with the iT-magnitude lim- 
its of SDF and SXDF. We find that the effect of limiting 
magnitude is not significant at Lho ^ 10"'^ '', 10^^'^, and 
2Q42.0 gj-g/g fQj. ^i^g redshift ranges of z = 0.5-1.0, 1.0- 
1.4, and 1.4-1.7, respectively, which are corresponding to 
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Fig. 6. The Ha luminosity functions of SDF galaxies (filled 
triangles) and SXDF galaxies (filled circles) in z = 0.5—1.0 
(top), z = 1.0-1.4 (middle), and z = 1.4-1.7 (bottom). 
The vertical lines indicate the luminosities corresponding to 
^Ha = 10~^^erg cm~^s~^ at each redshift interval. Estimates 
by previous studies are also plotted, whose symbols are indi- 
cated in the figure. The redshift ranges of the previous es- 
timates are as follows: Hopkins et al. (2000, HOO) in z = 
0.7-1.8, Yan et al. (1999, Y99) in 2 = 0.7-1.9, Tresse et al. 
(2002, T02) in z = 0.5-1.1, and Geach et al. (2008, G08) 
in 2; ~ 2.23, even though there is a small mismatch of red- 
shift range. The "high S/N "data of Hopkins et al. (2000) arc 
using only galaxies with spectroscopic confirmation or high 
S/N ratios. The data of G08 are obtained by a narrowband 
observation, while the others are spectroscopic observations. 
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Fig. 7. The A'-magnitude distributions of SDF (dotted) and 
SXDF galaxies (solid) in five Lhq (in ^rg/s) bins at 2 = 
0.5-1.0 (top), 1.0-1.4 (middle), and 1.4-1.7 (bottom). The 
vertical lines are the K limiting magnitudes of SDF (dashed 
line, So" at 2" aperture) and SXDF (solid line, 3cr at 2". 3 
aperture) . 



the Ha flux of Fhq > 10~^^erg cm~^s~^. These luminosi- 
ties are indicated by the vertical solid lines in Figure 6. 
Because galaxies for the BAO surveys considered in this 
work are brighter than this Ha flux, the incompleteness 
does not seriously affect the main results of this work. 
As another test, we estimate the number of Ha emission 
galaxies that are brighter than Fhq = lO^^^erg cm^^s^^ 
but have been missed in our sample because of the re- 
quirement of _fir-band detection, by using the i' selected 
sample. In this case, we cannot use the Bz'K diagram 
to remove stars and hence stellar contamination may in- 
crease, but we find that the number of bright emission 
galaxies increases at most 20% from our baseline analysis 
requiring if -band detection. 

The agreement with the previous data derived by spec- 
troscopic or narrow-band filter observations is reasonably 
good above the limiting magnitudes, providing some cre- 
dence to our study based on photometric line luminosity 
estimates. However, the brighter-end of the LF at z = 
1.4-1.7 estimated by the SDF data is significantly lower 
than those by the SXDF data and previous studies, which 
cannot be explained solely by the limiting magnitude ef- 
fect. We find that the reason for this discrepancy is con- 
fusion between galaxies at z ~ and z 1.4 due to the 
lack of the mid-infrared data in 3.6 and 4.5 /im bands in 
the SDF. We confirm this by checking that the LFs esti- 
mated by SXDF galaxies without using mid-infrared data 
are similar to those by SDF galaxies, thus the SXDF LF 
is preferred in the bottom panel of Figure 6. 

^.2. [Oil] Luminosity function 

We also calculate the [Oll]A3727 luminosity function 
in SDF and SXDF, which are shown in Fig. 8, as 
well as those by previous studies (see figure caption). 
The limiting magnitude effect is analyzed in the same 
way as with the Ha luminosity function, by breaking it 
down in if-magnitude bins (Fig. 9). The figure shows 
that the limiting magnitude effect is not significant for 
L[oii] ~ 10^^ ■^ 10^^ ■^ and lO^i-^ erg/s at z= 0.5-1.0, 1.0- 
1.4, and 1.4-1.7, respectively, which are corresponding to 



[Oil] 



> 10 



-16.5gj.g 2g 



These luminosities are also 



shown by the vertical lines in Figure 8. As in the analysis 
of the Ha LF, we estimate the number of [Oil] emission 
galaxies brighter than f[on] ~ 10^^^'^erg cm~^s~^ using 
the sample selected by i-band detection regardless of the 
i^T-band detection, and we find that the increase of emis- 
sion galaxies is at most ^ 30%. 

Although the results from different studies are reason- 
ably consistent at z = 0.5-1.0 and 1.0-1.4, the agreement 
is rather worse than for the Ha LF. The reason could be 
due to the large scatter of the L[on]^SFR relation or un- 
certainties in our calibration procedure. The paucity of 
galaxies at the faint end of our LF at z = 1.0-1.4 com- 
pared with previous studies is most likely due to the limit- 
ing magnitude effect. However, an even larger discrepancy 
is found at the faint-end of the LF at z = 1.4-1.7, which 
can not be explained solely by the limiting magnitude ef- 
fect. (The disagreement between the [Oil] LFs of SDF 
and SXDF is again caused by the lack of the mid-infrared 
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Fig. 8. The same as Fig. 6, but for [On] line luminosity 
functions. The vertical lines indicate the luminosities corre- 



sponding to Fi 



[Oil] 



IQ- 



-I6.5gj.g pj^— 2g-l each redshift 
interval. The redshift ranges of the data of previous studies 
are as follows: Hippelein et al. (2003, H03) at 2 ~ 0.88 (top 
panel) and ~ 1.18 (middle); Ly et al. (2007, L07) at z ~ 0.89 
(top), 2 ~ 1.19 (middle), and z ~ 1.47 (bottom); Takahashi et 
al. (2007, T07) at 2 ~ 1.19 (middle); Zhu et al. (2008, Z08) 
at 2 = 0.752-0.926 (top), 2 = 1.108-1.277 (middle), and 2 = 
1.277-1.446 (bottom). The data of H03 and Z08 are obtained 
by spectroscopic observations, while those of L07 and T07 are 
based on narrowband observations. 



Fig. 9. The same as Fig. 7, but for [On] line luminosity 
functions. 



No. ] 

data in SDF.) One possibility is the contamination by dif- 
ferent redshift emission line objects in the sample of Ly 
et al. (2007) which was selected by narrow-band photom- 
etry. However according to Ly et al. the contamination 
rate of their [O II] sample, estimated from their spectro- 
scopic sample, is only 10-20%. The cosmic variance due 
to the small survey volume of Ly et al. may also have 
caused the discrepancy. Further systematic uncertainties 
in our photometric line luminosity estimates, such as the 
photometric redshift errors at z > 1.4 (see Fig. 1), may 
also be the origin of the discrepancy. Since the Ly ct al. 
data is the only previous study that can be compared with 
the faint-end of our LF in this redshift range more data 
from independent studies are highly desirable. 

5. Properties of Emission Line Galaxies 

In this section we discuss implications for future BAO 
surveys, based on the emission line luminosity estimates 
obtained above. We intend to make our results as gen- 
eral as possible, but a specific discussion assuming exist- 
ing proposals for future BAO surveys is also useful and 
illustrative. In particular, we discuss the proposed BAO 
surveys using FMOS and WFMOS as examples of Ha and 
[On] surveys, respectively. We consider three representa- 
tive redshift ranges in the following analysis: 0.5 < z < 1.0, 
1.0 < z < 1.4, and 1.4 < z < 1.7, as used in the LF esti- 
mates. The whole redshift range 0.5 < z < 1.7 roughly 
corresponds to the range that can be probed by optical 
[O II] and near-infrared Ha surveys. 

5.1. Galaxy Surface Density as Functions of Line Flux 

First we check the available number of emission line 
galaxies in the three redshift intervals. These are the max- 
imum numbers of emission line galaxies that can be used 
in any BAO survey. We plotted the cumulative number 
of SXDF galaxies as a function of Ha and [Oil] flux in 
Figure 10. The solid horizontal line in the panel indicates 
the surface number density of nP = 1 assuming z ^ 1, 
where n is galaxy number density and P is the ampli- 
tude of the power spectrum of galaxy density fluctuation 
at fc = 0.20/iMpc-\ which is roughly the linear-nonlinear 
transition scale at z 1 . The survey power is maximized 
by choosing nP ~ 1 when the observing time is fixed (Seo 
& Eisenstein 2003; Glazebrook & Blake 2005). Here we 
used the formula of Eisenstein & Hu (1998) to calculate 
the CDM power spectrum. The linear clustering biases 
b = 1.19 and 1.17 are taken into account (Pgai = b^PcBu) 
for Ha and [O II] lines, respectively. These values are from 
the bias estimates for emission line galaxies at 1.0<z<1.4 
by stellar mass estimates, as described in §5.4. 

We define i%Q,min as the flux above which the number 
of galaxies is equal to the surface number density cor- 
responding to nP = 1, and hence -FHa.min gives the mini- 
mally required flux sensitivity to perform any BAO survey 
by emission line galaxies. We find that these values are 
log(FHa,min/[crg cm~^s~^]) = —15.5, —15.5 and —15.8 for 
the redshift intervals of 0.5 < z < 1.0, 1.0 < z < 1.4 and 
1.4 < z < 1.7, respectively. These fluxes are brighter than 
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Fig. 10. The cumulative number counts of galaxies as func- 
tions of Ha (top) and [O ll] (bottom) flux. The solid, dashed, 
and dotted lines are for the redshift ranges of 0.5 < z < 1.0, 
1.0 < 2: < 1.4, and 1.4<z<1.7, respectively. The solid horizon- 
tal lines indicate the surface number density corresponding to 
nP = 1 at z ^ 1. (This density does not change significantly 
for the three redshift ranges.) The dashed horizontal lines 
represent the FMOS and WFMOS fiber density in the top 
and bottom panels, respectively. 

the FMOS sensitivity, ~ log(FHa/[crg cm'^s-^]) = -16.0 
for 1 hour integration (Eto et al. 2004), and hence there 
are a enough number of Ha-emitting galaxies for a BAO 
survey by FMOS. The dotted horizontal line represents 
the fiber density of FMOS (400 in 0.2 deg^), which is close 
to the galaxy number density corresponding to nP = 1, 
meaning that FMOS is an efficient instrument for a BAO 
survey than can collect a necessary number of galaxy spec- 
tra by one time observation per field, although BAO was 
not the original scientific target of this instrument. 

Similarly, we define i^[oii],min for [On] line, and we 
find that these are log(F[oii],min/[crg cm~^s~^]) ~ —16.1, 
— 16.1 and —16.4 for the three redshift intervals, respec- 
tively, which should be compared with the line sensitivity 

of WFMOS, log(F[oii]/[erg cm-2s-i]) 16.5, for 0.5 

hour integration (Parkinson et al. 2007). The horizontal 
dotted hue is the planned fiber density of WFMOS (3000 
in 1.8 deg^). As in the case of Ha line for FMOS, there is 
a sufficient number of emission line galaxies that can be 
used for a [Oil] BAO survey by WFMOS. 
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Fig. 11. The cumulative galaxy numbers of the Ha emitting galaxies in SXDF as a function of magnitudes in popular band filters. 
They are shown separately for the three redshift intervals as indicated. The dashed, solid, and dotted lines show the number of 
emission line galaxies brighter than the three different threshold Ha fluxes of 10~^^'^, 10~^^, and 10~^®'^erg cm~^s~^, respectively. 
The horizontal lines (solid and dashed) indicate the surface number density corresponding to nP = 1 at each redshift interval and 
the FMOS fiber density, respectively. 
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Finally we compare this result with that of the WiggleZ 
survey, a BAO survey targeting emission-line galaxies in 
the rcdshift range of 0.2 < z < 1.0 (Glazebrook et al. 
2007; Blake ct al. 2009). The number of galaxies brighter 
than -F[oii] = lO^^^erg cm^^s"^ at z =0.6-0.7 by our es- 
timates is 210 deg~^, while that of actually observed by 
the WiggleZ survey is 35 deg~^. The paucity of galaxy 
number of the WiggleZ is likely due to strong UV and 
r-band selection criteria of the survey. 

5.2. Galaxy Surface Density as Functions of Broad-band 
Filter Magnitudes 

Even if there are a sufficiently large number of galaxies 
for BAO surveys, they must be pre-selected by imaging 
surveys, and the required depth for such photometric se- 
lection is crucial in planning BAO surveys. Figure 11 
shows the cumulative galaxy number of the Ha emitting 
galaxies in SXDF function of magnitudes in popu- 
lar band filters. They are shown separately for the three 
redshift intervals and different Ha threshold fluxes as in- 
dicated. Figure 12 is the same as Figure 11, but for [Oil] 
and WFMOS. 

From these figures, one can infer the necessary (but not 
sufficient) depth in each band filter for photometric selec- 
tion of targets. For example, we need imaging survey lim- 
its deeper than B = 23.9, i' = 23.3, z' = 22.8, J = 22.3, and 
i^T = 21.7 to detect BAO target galaxies brighter than the 
Ha sensitivity of FMOS (> 10-^^ ° erg cm'^s'^) with a 
number density of FMOS fibers at 1.0 < z < 1.4. Similarly, 
the required imaging depths for the WFMOS [O II] sensi- 
tivity (> 10~i6-^ erg cm-^s-i) become B = 23.8, i' = 23.3, 
z' = 22.9, J = 22.4, and K = 21.9. The simflar magnitudes 
for FMOS and WFMOS indicate that the BAO targets for 
these two surveys have similar SFRs. While Ha line lu- 
minosity is generally larger than [Oil], NIR observations 
are less sensitive than optical observations, and hence the 
required integration time for Subaru becomes comparable 
for Ha observed by FMOS and [O II] observed by WFMOS 
for a given SFR. 

5.S. Stellar Mass, Extinction, and Size 

Here we discuss some physical properties of the emission 
line galaxies that are useful for a planning of a BAO sur- 
vey. First, we show the stellar mass distribution of SXDF 
galaxies that have been calculated in the photometric rcd- 
shift estimates, in Figure 13 with the threshold fluxes of 
10~^^ erg cm~^s~^, for Ha, or 10~^^'^ erg cm~^s~^, for 
[On] flux, in the three redshift intervals. We compare 
the mass distributions of these results with those of all 
galaxies in SXDF, and then we flnd that almost all galax- 
ies which have small masses are emission line galaxies. 
The logarithmic means of the stellar mass distributions 
are presented in Table 2, for several different line flux 
thresholds. The mean stellar mass does not change sig- 
nificantly with the threshold line flux in a fixed redshift 
interval, and does not change either with redshift for a 
fixed threshold flux. The mean stellar mass of line emit- 
ting galaxies is rather low, but we note that there is quite 
a large spread of stellar masses from very massive galax- 
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Fig. 13. The stellar mass distributions of Ha and [0 11] omis- 
sion galaxies in SXDF at z =0.5-1.0 (top), 1.0-1.4 (middle), 
and 1.4—1.7 (bottom). The Ho and [On] threshold line fluxes 
are 10~^^ and 10~^^'^ erg cm~^s~^, respectively. We also 
show the mass distributions of all galaxies in the SXDF sam- 
ple in the same redshift range, for the comparison with the 
results of emission line galaxies. 

ies down to dwarfs. This result indicates that the line 
luminosity is not strongly correlated with stellar mass of 
galaxies. 

Next we examine the extinction (Ay) distribution of 
SXDF galaxies estimated by the photometric redshift 
method in the three redshift ranges, as shown in Figure 
14. The extinction is distributed with a large scatter, with 
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Table 2. Logarithmic means of the stellar mass of SXDF galaxies 



Line flux [ erg cm ^s' 



0.5<z< 1.0 



logAf, [ Mo ]* 
1.0 < z < 1.4 



lA<z< 1.7 



log^^Ha > -16.5 
logi^Ha > -16.0 

logJ^Ha > -15.5 



9.8 
9.9 
9.7 



9.7 
9.6 

9.4 



9.9 
9.6 
9.6 



logF[o„] > -17.0 
logFfoii] > -16.5 



9.6 
9.7 
9.5 



9.6 
9.5 

9.3 



logfronl > -16.0 
' yec J^'ig. 13 tor the presentation of stellar mass distributions tor threshold 
fluxes of logFHc > —16.0 and logFjon] > —16.5. 



9.8 

9.5 
9.7 



mean values of Ay — 1.2 and emission line galaxies have 
larger extinctions in this figure. However the distribution 
of galaxies selected by [Oil] at the lowest redshift range 
of 0.5 < z < 1.0 shows a significant number of galaxies 
with Ay ~ compared with other redshift ranges. It is 
uncertain whether this is real or caused by systematic un- 
certainties of photo-z calculation. 

We also calculate the histogram of FWHM sizes of the 
SXDF galaxies in the three redshift ranges, as shown in 
Figure 15. One can compare these sizes with the fiber 
aperture or slit width in a planned survey to estimate 
the loss of emission line light due to the extended size of 
galaxies. The FWHM size was measured in the B band, 
because the galaxy profile in the shortest wavelength is 
expected to be close to that in emission lines. Seeing 
0.8 arcsec FWHM) has not been subtracted, mean- 
ing the values correspond to the actual image size under 
realistic observing conditions. In Figure 15, every distri- 
bution of FWHM size has a lower cut-off at approximately 
0.75 arcsec and gradually decays with increasing size. The 
lower cut-off is due to the original image quality of the op- 
tical images in SXDF. For reference, the fiber diameter of 
FMOS is 1".2 diameter, and the planned fiber diameter 
of WFMOS is also about 1". These results indicate that a 
significant fraction of emission line flux could be lost out 
of the fibers in large size galaxies, and this effect must be 
taken into account in BAO survey planning. 

5.4- Linear Bias 

The bias of clustering is important to estimate the de- 
tectability of BAO signatures. We therefore calculate lin- 
ear bias of clustering by using the baryon mass estimate of 
galaxies in SXDF. The baryon mass is estimated from the 
photometric redshift calculations, including stellar mass 
and gas mass that is not yet converted into stars. (In the 
case of the SED template of constant SFR, we suppose 
that gas mass is equal to the stellar mass estimated from 
the phot— z calculations.) The estimated mean gas mass 
is consistent with that of UV-selected galaxies at z ~ 2 
obtained by Erb et al. (2006). We then make a rough 
estimate of dark matter halo mass, assuming the cosmic 
baryon-to-DM ratio, Jlb/riDM: 

AfDM = ^(A/,+Mga.), 
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Fig. 14. The same as Fig. 13, but for the extinction (Ay) 
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Fig. 15. The same as Fig. 13, but for the B-band FHWM 
size distributions of Ha and [O ll] emission galaxies. The see- 
ing is not deconvolved in these sizes. 
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Fig. 16. The top panel shows the linear bias estimated for 
the SXDF galaxies that are brighter than a given Ha flux, 
for the three redshift intervals as indicated in the panel. The 
bottom panel is the same as the top, but for the [O ll] line. 

where fluM = — ^b, and and fib are the density pa- 
rameters for matter and baryons, respectively. We then 
calculate the logarithmic mean of dark halo mass for a 
given threshold line flux in a given redshift interval, and 
then calculate linear bias from the estimated mass us- 
ing the bias model of dark matter haloes of Cooray & 
Shcth (2002). The biases estimated in this way are shown 
in Figure 16 for the SXDF galaxies as functions of the 
threshold Ha and [O II] fluxes. The estimated linear bias 
is typically b ~1-1.5. Since the stellar mass of galaxies 
does not change with the threshold flux or redshift, the 
linear bias does not depend strongly on the threshold flux, 
but it becomes larger at higher redshifts. 

We estimated the linear bias assuming that the DM- 
to-baryon mass ratio of galaxies is the same as the cos- 
mic value, but it may not be the case, and it is likely an 
underestimation because of the following reasons. First, 
the gas mass estimated by the star formation history of 
photo-z templates (exponential SFR evolution) is highly 
uncertain, and gas mass can be much larger than stellar 
mass in young, line-emitting starburst galaxies (Erb et al. 
2006). Second, we have implicitly assumed that all the 
baryons are either in stars or in cold gas participating in 
star-formation. However, a significant amount of diffuse 
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hot gas could also lie in dark halos. It has been found 
empirically that the ratio of dark halo and stellar mass is 
approximately 2-20rJDM/^^ti for different luminosities of 
galaxies at 2; ~ 2 (e.g., Hayashi et al. 2007). Therefore we 
also calculate the biases with a 3 and 10 times larger DM 
mass than the above estimates, but the linear biases are 
only modestly increased by 10% and 30%, respectively. 

We can compare the bias estimates with those by 
clustering analysis in the literature at z ~1. 0-2.0, by 
converting correlation length into bias of line-emission 
galaxies. Blake et al. (2009) obtained h 1.6 for [Oil] 
emitters at z =0.5-1.0 from the real space correlation 
function of the WiggleZ project data (Glazebrook et al. 
2007; Blake et al. 2009), with a threshold [Oil] flux of 1 
hour integration on the 3.9m Anglo- Australian Telescope 
(~ 10~^^erg cm~^s~^). This bias is larger than our esti- 
mate, and this can not be explained by increasing the DM- 
to-baryon mass ratio of galaxies. However, this is likely to 
be due to the strong additional r-band selection criteria 
of the WiggleZ survey (20 < r < 22.5). We found that, 
when we add 20 < r < 22.5, our bias estimate is increased 
by 20% as well. Therefore, our bias estimate is in rough 
agreement with the WiggleZ clustering data when the in- 
crease of our bias prediction by larger SIdm/^^!) values is 
taken into account. Geach et al. (2008) obtained h = 2.0 
for Ha emitters at z '^2.2 by angular correlation function, 
with a threshold Ha flux of ^ 10""'^^ erg s~^ cm~^. This is 
out of the redshift range that we investigated, but it seems 
also roughly consistent with our estimate considering the 
increasing trend of bias with redshift. 

5.5. Clustering Property 

Next wc examine the ficld-to-field variation of the ob- 
served number of emission-line galaxies in SXDF, and 
compare it with those expected from the linear bias es- 
timated above and the CDM structure formation theory. 
A detailed study on the clustering properties by using an- 
gular correlation function will be reported elsewhere. We 
divided the region in SXDF used in this work (the region 
having both optical and NIR data) into four sub-regions 
(SXDF 1-4) with roughly the same area of 0.18 deg^ 
(see Figure 17 and Table 3). 

We then estimate the standard deviation cTgai.obs of the 
galaxy density fluctuation ^gai = Spgai / Pgai on the volume 
scale of the sub-fields as follows. When the surface galaxy 
number density n,; in the i-th sub-field obeys a distribution 
with the mean value n, the unbiased estimate for the stan- 
dard deviation (T„ of the distribution from the observed 
number densities is given by: 



N 



1 ^ 



where nobs = /N is the observational estimate of 
the mean value. Then, (Tgai,obs is simply estimated by 
CTn/fiobs- The observational error in the 

^gal,obs estl 

mate can be calculated by approximating that '^{ni ~ 
nohsY /^n obeys to distribution with — 1 degrees of 
freedom. The value of CTgai,obs and its error are thus esti- 



SXDF 1 



SXDF 3 



SXDF 



SXDF 



Fig. 17. Overview of SXDF. The regions indicated by thin 
solid lines are the five fields of view of Subaru/Suprime-Cam 
(SXDS-C, N, S, W, and E). The region shown by dashed lines 
is the area where NIR data by UKIRT/WFCAM is available. 
All of the field is covered by mid-infrared observations by 
Spitzer/IRAC. We used the overlapping region of the optical 
and NIR data, and the thick solid lines show the four sub- 
fields defined in this work. 



mated for galaxies brighter than 10~"'^^erg cm~^s~'^ (Ho;) 
and 10~^^ '^erg cm~^s~^ ([Oil]), which are given in Table 
3 in the three redshift ranges. 

To compare with these results, we calculate the theoreti- 
cally expected standard deviations, 6ctcdm, where ctcdm is 
the standard deviation of dark matter fluctuation and b is 
the linear bias of line emitting galaxies. The value of (Tcdm 

is calculated assuming a spherical top-hat window func- 

1 

tion with radius R, and R is assumed to be i? = (31//47r) ^ , 
where V is the averaged volume of the four sub-flelds of 
SXDF. The results are given in Table 3, to be compared 
with (Tgai,obs. The values of 6(Tcdm are in agreement with 
Cgai.obs in 1.0 < z < 1.4 and 1.4 < z < 1.7, but in disagree- 
ment in 0.5 < z < 1.0. The origin of the discrepancy only 
in 0.5 < z < 1.0 is not clear, but the estimate of (Tgai,obs 
is based only on four subflelds. We flnd that the chance 
probability of getting this level of deviation is about 4% 
and 20% for Ha and [O II] line cases, respectively, and 
hence the possibility of a statistical fluke cannot be ex- 
cluded. 



(1) 6. Example of Color Selection 



Actual selection procedures in a particular BAO sur- 
vey will depend on various conditions that are unique to 
the survey, such as available band fllters of input imag- 
ing surveys. Therefore it is difficult to derive generally 
useful results, but here we test some simple two-color se- 
lection methods to select emission line galaxies brighter 
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Table 3. Field-to-field variation of galaxy number 



Field 


Area [deg^] 






Galaxy Number 


in Each Field 










lOf 

z=0. 5-1.0 


2=1.0-1.4 


16.0 

2=1.4-1.7 


logi^[oii]*> -16.5 

2=0.5-1.0 2=1.0-1.4 2 


= 1.4-1.7 




1 7Q 


i- i-Ou 




oou 


1523 


1238 




SXDF 2 


0.199 


1235 


1311 


371 


1768 


1241 


322 


SXDF 3 


0.176 


1119 


1044 


419 


1610 


1005 


358 


SXDF 4 


0.178 


1134 


907 


394 


1691 


850 


328 


'^gal,obs 




0.012 


0.168 


0.099 


0.047 


0.153 


0.093 


^f^gal,obs ^ 




±0.009 


±0.137 


±0.081 


±0.038 


±0.124 


±0.075 






0.125 

-ii 


0.105 


0.104 


0.120 


0.103 


0.109 



* In units of line flux [erg cm 

t The standard deviation of galaxy density fluctuation, 5pgai/Pgal estimated from the observed galaxy numbers in sub-fields. 



-f The observational error of (Tgai.obs- 

^ The expected density fluctuations predicted from the CDM theory and the bias estimates, which should be compared with 



than Fhq > 10~^^erg cm~^s~^: one is using only optical 
photometries {Bi'z'), and the other is including the NIR 
K band {Bz'K). We do not show in detail the results for 
[O II] , but we have confirmed that the results arc similar 
when one considers the corresponding [O II] line flux of 

^Q-16.5 gj.g cn2-2g-l^ 

Before we examine the two-color diagrams, we should 
select relatively bright galaxies as potential targets for 
BAG surveys, using photometric information. From the 
magnitude distribution of emission-line galaxies (Figure 
11), we set the magnitude threshold as S < 24. We then 
show the color-color diagrams of all the SXDF galaxies 
with i? < 24 for the two cases, in Figure 18 and Figure 19. 
Galaxies that arc brighter than Fhq = 10~^^erg cm~^s~^ 
are shown by colored dots, for the three different redshift 
intervals. From these results, we define the color selection 
criteria to efficiently select the emission-line galaxies, as 
shown by the box in the same colors for the three red- 
shift intervals. These criteria were chosen not only by 
the number of target emission line galaxies, but also by 
the expected success rate examining the contamination of 
non-emission line galaxies (black dots) . These criteria can 
also be applied for [Gll] emitters brighter than 10~^^'^ 
erg cm^^s"^ because the distributions on the two-color 
planes are similar to those in the Ha cases. 

Figures 20 and 21 show the number of emission- line 
galaxies brighter than three threshold Ha line fluxes of 
]^g-i5.5^ 10~^^, and 10~^^-^crg cm~^s~\ as a function of 
the threshold B magnitude when we select galaxies by 
these color criteria. It can be seen that for the two low 
redshift bins sufficient galaxies {nP ^ 1) are obtained at a 
depth of B = 24 with a 10~^^erg cm~^s~^ threshold. For 
the highest redshift bin we need to go deeper {B = 25) to 
attain sufficient galaxies. 

The success rate of the selection, i.e., the fraction of 
galaxies having larger emission line flux than the thresh- 
old and in the correct redshift range is also shown in these 
figure, compared with the total number of all galaxies sat- 
isfying the B magnitude threshold and color selection cri- 
teria. These results indicate that typically a success rate 
of about 40% in the two low redshift bins can be achieved 



1 




1.4<z< 1.7 

-0.2 ' ^ ' ' ' 1 

0.5 1 1.5 2 2.5 3 

B - i' [mag] 

Fig. 18. The Bi'z' plot for SXDF galaxies with B < 24. 
Galaxies with Ha flux brighter than 10~^^ erg cm^'^s"^ in 
the three redshift ranges are indicated by red, green, and blue 
dots, while the other galaxies in black dots. The color selec- 
tion criteria used for the three redshift ranges are indicated 
by the corresponding colors. 




0.5 1 1.5 2 2.5 3 



B - z' [mag] 

Fig. 19. The same as Figure 18, but for Bz'K diagram. 
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[Vol. 



0.5<z< 1.0 



1.0<z< 1.4 



1.4<z< 1.7 




Bab [mag] 



Fig. 20. The counts and probability of successful selection in the three redshift intervals by the Bi' z' selection method, for galaxies 
brighter than a given threshold B magnitude (abscissa). The galaxy sample is SXDF. The horizontal lines (solid and dashed) in 
the top panels indicate the surface number density corresponding to nP = 1 at each redshift interval and the FMOS fiber density, 
respectively, 
respectively. 
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Fig. 21. The same as Fig. 20, but for the Bz' K selection. 
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Fig. 22. The same as Fig. 20, but for the Bi' z' selection on the 1.0 < z < 1.7 joint condition. The left, middle, and right panels 
show the counts and probability of the successful selection for galaxies in the redshift intervals of 1.0 < z < 1.4, 1.4 < z < 1.7, and in 
the joint redshift interval of 1.0 < 2 < 1.7. 




Fig. 23. The same as Fig. 22, but for the Bz' K selection on the 1.0 < z < 1.7 joint condition. 
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by this simple color selection for the threshold line flux of 
Fna = lO^^^erg cm^^s^^. For the case for [Oil], we also 
find that typically a success rate of about 40% in the two 
low rcdshift bins can be achieved for the threshold line 
flux of i^[oii] = lO^iS-^erg cm-^g-i. 

Finally we consider the two-color selection for high red- 
shift galaxies at z = 1-1.7 using Bi'z' and Bz'K pho- 
tometries. This may be effective if we do not discriminate 
galaxies beyond or below z = \A and are simply interested 
in selecting galaxies at z > 1. We define the selection crite- 
ria for Bi' z' and Bz'K as shown by the two high redshift 
boxes in Figure 18 and 19, respectively. From Figure 22 
and 23, we can select a sufficient number of galaxies at 
1.0 < z < 1.7 by a survey deeper than B = 24 with a success 
rate of 50-60%. If the survey depth is deeper than B = 25, 
we can select a sufficient number of galaxies both in the 
two redshift intervals of 1.0 < z < 1.4 and 1.4 < z < 1.7. 

We find that both Bi'z' and Bz'K selection deliver 
comparable performance for selection in terms of mag- 
nitude depth required and success rate. The choice would 
then come down to the ease of obtaining the observational 
datasets. 

Possible improvements would be to try three-color 
selection boxes, or multi-parameter methods in an N- 
dimensional photometric space. These improvements are 
beyond the scope of this paper but we will explore such 
possibilities in future work. 

7. Summary 

We made an estimate of Ha and [Oil] emission line 
luminosities of SDF and SXDF galaxies in a total field 
area of 0.846 deg^ by a photometric redshift based ap- 
proach, for the purpose of studying photometric selection 
procedures of high redshift emission line galaxies useful 
for future spectroscopic BAO surveys. Our line luminos- 
ity estimates were tested and calibrated by using the SDSS 
galaxies with spectroscopic line luminosity measurements 
at low redshifts. The line luminosity functions of Ha and 
[On] were thus derived in three reshift intervals in z = 
0.5-1.7, and they are in reasonable agreement with the 
previous studies based on spectroscopic or narrow-band 
filter observations. 

We examined the number density and properties of 
emission line galaxies that can be used for future BAO 
surveys, taking FMOS and WFMOS instruments for the 
Subaru Telescope as representative examples. We have 
confirmed that there are sufficient number of emission line 
galaxies to perform BAO surveys in the redshift ranges 
0.5 < z < 1 and 1 < z < 1.4 which are bright enough to 
be detected by 8-10 m class telescopes with exposures 
shorter than one hour. We showed cumulative distribu- 
tions of magnitudes in popular band filters for emission 
line galaxies with line fiux brighter than given thresholds 
in a range of redshift intervals, and these can be used to 
estimate the required depth of an imaging survey from 
which spectroscopic targets are selected. 

We also estimated the stellar mass distributions, extinc- 
tion distributions, FWHM size distributions, and linear 



biases of line emitting galaxies. The linear biases were in- 
ferred from the stellar mass and redshift estimates, com- 
bined with the standard theory of structure formation. 
The inferred biases are consistent with the field-to-field 
variation of galaxy numbers in four subfields in SXDF, 
and also with the estimates based on clustering analysis 
for similar galaxies in previous studies. 

Finally, we showed some examples of color selection of 
target galaxies for spectroscopic BAO surveys, based on 
simple color-color selections {Bi'z' and Bz'K). It is found 
that target galaxies for BAO surveys can be selected with 
a reasonable accuracy and success rate. We find that 
with this selection for z < 1.4 we need to reach a typi- 
cal depth of _B = 24 (with other filters being required for 
color-selection) and for 1.4 < z < 1.7 we need to reach 
B = 25. We find the Bi'z' and BzK' selection methods 
give comparable efficiency. In general we find that it is not 
possible to efficiently select galaxies in the high-redshift 
1.4 < z < 1.7 using only two-colour selection. 

These results give guidance for the planning of future 
BAO surveys and demonstrate that 8-m telescopes with 
spectroscopic fields of view of order one degree (instru- 
ments such as FMOS and WFMOS) would be effective 
machines for carrying out such large cosmological surveys. 

The data of photometric magnitudes, redshifts, and line 
luminosities used in this work are available on request to 
the authors. 
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